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INTRODUCTION 

The method of characteristics for three dimensional axially 
symmetric bodies was used to determine the velocity distribution 
about the nose of the Corporal E rocket, a rocket projectile, for 
fdach numbers 2 and 4. From the velocities the pressure distrib- 
ution was determined ola a drag coefficient computed. 

Por a starting point the nose of the vrojectile was apvroxe- 
imated for a short distance by a cone and the Taylor-laccoll"” 
solution to this problem was used. This solution gave the anzle 
of shockwave and the body. The Sauer graphicelenumerical itera- 
tion method was used for the remainder of the golution. 

Yreliminary calculations and the work for /iach number 3 were 


carried out by Dr. H. K. Forster to whon we aro indebted for 


instruction and aid with this work, 


“Taylor, @. I. and Maccoll, J. W., Proc. Roy. Soc., A, VYOl. 139, 
pn 278-311 1933. 





CHARACTERISTICS® 
The following discussion is for two dimensional flow and igs 
taken from "Introduction to Aerodynamics of a Compressible Fluid", 
by H. #. Liepmann and Allen E. Puckett. 
The term "characteristic" is related by definition to proper- 
ties of certain partial differential equations. Yor any second 


order, partial differential equation of the tyne 





a x2 2 
(1) A Se t 6B ig, +O yy 


a2 yop hz 
ee ee 
two families of curves may be defined by the total differential 


equations 


(2) oe, . Bae = Ac 


d x A» 





where A, B, and C are functions of x and y. These curves (2) are 
known as the characteristics of equation (1). 

These characteristics heve certain proverties which make then 
useful in an approximate solution process. These properties will 
be listed without proof: 

1) They are inveriant under transformation. 

2) A solution Z(x,y) to equation (1) may be composed of 
separate integral surfaces joining continuously with continuous 
derivatives along a line. This line is a characteristic. This 
property of the branching of integral surfaces along a character-~ 
istic permits us to regarc the line as a wave front along wnich a 
disturbance is propageted through a fieclc. 


*Lieomann, H. W. and Puckett, Allen ©. Introduction to Aerodynamics 
of a Compressible Fluid, wn 2292232 1947 


Se 
3) Boundary conditions alone any line, C, in the field will 
define 2 golution within a triangle bounded by © and the cherac- 
teristics through the end points of C, Boundary conditions along 
a characteristic, however, will not define a solution in any region; 
boundary conditions must be given along two intersecting character- 
istics of different families. The exact equation of the velocity 


potential is: 
: 2 2% 
(3) (a* -u)2£ 7 pans f f + (a2-y) 2 =O 
J 


yy 
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dX na dy 
If a new dependent variable is defined 
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x» 
P 


where AL = 


then when a transformation is made to the hodograph plane, an ele- 
ment of the F surface in the (u,v) plane is defined by an elenent 
of the (x,y) surface. (Legendre Transformation). The transform 


ation of equation (3) becomes 


(4) (a*-v*)E, tua, + (@-v7)F, =O 


The characteristics of the equation are 


Oat POEL ED 


(5) dv a as pt ar? (a *~ra*}(a*-av*) 
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This equation can be integrated immediately, "a" having been expressed 
as a function of u and v. Thus in the two dimensional case a charac- 
teristic chart may be constructed and is known as the "Prandtl- 
Busemann Characteristic Diagran." 

When an attempt is made to use this method in three dimensional 


problems it is found that it cannot be done. A Legendre Transformation 


a. 


leads to a non-linear equation in this case. In general, in three 
dimensional flow, it is not nossible to construct a char<cteristic 


diagram. Insterd, a senarate construction must be carried out for 


each floy, 


SAUER MBPHOD 
The 3auer solution of supersonic axially symmetric floyv is a 
erenhicelenumerical iteration method. 
"or a non-viscous, stendy, isentropic flow of a verfect fluid 


without body forcds, tne equation of notion mey be obtained by 





eliminating the pressure end density from the force equations end 
the equation of continuity. Transformed to cylindricel coordinates, 


end assuming axial symmetry, this becomes 


. +) a £ 
 G- AE + YE - a i the +E 28 


Assuming irrotational floy the velocity commonents may be re vlaced 
by the proper derivatives of the velocity potential gp . 

In the axially symmetric case, given a velocity vector, VY, one > 

rthogonal 

can introduce in en axial nlane a set offaxes in which X!' is elong 
the velocity vector. Then u' = V and v'! = O where u! is the velocity 
in the X' direction end v! is the velocity in the Y' direction. 

If at a particular noint the X%', Y' coordinates are introduced 
equation (6) becomes 

; / 
_ a Be (oe ~ = a) 

and introducing j » tor irrotational flow 


(7) (i- 4% ag , 1g jae a, 


ai oxi” J ae A i” ae 


~S 


The Mach angle is defined as oo: a! = where a is the local 
sneed of sound and V the local fluid velocity. <A Mach line is a 
line inclined to the direction of fluid flow by an angle equal to 


the ilach angle. 


At a point, where the velocity is known tyvo Mach lines may 
ve drawn. In this case vhere the velocity is along the X!' axis, 
one Mech line will ve at an engle and one et an angle fron 


theme «epmies §6€6(See Pie. 1 a). 





Fig. (1) a 


These tach lines, which may be identified as the character- 
istics of eouetion (6) define a new set of coordinate axes, and 
From the definition of the Mech anrle, sin it is seen that 


substituting this in equation (7) we ret 


oe 
TON 


(8) cot *« a 


© pb 
see a = 





Fig. (1) d 
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sudstituting these relations in (8) 
ee 


wep AOA 
Phis equation can be written 

J (4 > . ee 
(9) is 7 

4, ( 2) = a J Aon A ae 


where p oa q indicate c projections of the velocity vector along 


the } and n directions respectively. fo 


fenlacing the derivatives of (9) by a quotient of increments 
and clearing fractions the equations become 
dg = SS ee df 
(10) = 
dp - YY new x dN 
a 
being the inerements in the directiona of the Mach lines 


and 
and dp and dq the projections of the velocity increment unon the 


Mach lines. 





PROCKDURE 

The nroperties of the characteristica enable us to nake an 
annroximate solution. The mechanical vrocess is a simultaneous 
construction of velocities in the hodogreph rlene and of /l‘ach 
lines in the nhysical plane. Step by sten smell areas in the 
pnysical nlane correspond to points in the hodograph plane since 
the flow field in the physical nlane is divided into a finite 
number of small regions bounded by Mach lines over hich averages 
are taken for comruting velocities in the hodograph vlane. Hence 
the accuracy of the computation will denend on the size of the 
regions bounded by Mach lines. The boundary conditions in the 
ohysical plane determine the direction of tho velocity in the 
hodogravh plane; at the shock wave the magnitude is aleo deter- 
mined. Tha points in the hodograph plana when found by an iter- 
ation process based on estimetes will in turn determine the Mach 
lines in the pnhysicel vlane. 

The two boundaries which define the boundary conditions are 
the shock wave attached to the nose and the surface of the body. 
Between these boundaries a net of Hach lines is constructed by 
proceeding from point to point until, in the entire ares about 
the body and at voints along the surface, the velocities are 
determined. This Mach net is shown in Fig. (5), Fig. (6) or Fig. 
(2). 

The angle of the shock wave -t the nose is obtained fron 
Taylor-Meccoll solution of flow pee) ® cone with the same nose 
angle. For this purpose, denending upon the geometry of the body 


the assumption that the nose is conical will be good at least a 


a. 
small distance back from the vertex. Taylor and Maccoll showed 
that velocity megnitudes and directions ere constant in super- 
sonic flow about a cone alone radial lines through the vertex. 
This solution furnishes as many -oints 2s needed to start the 
cheracteristic method solution. These points can be taken close 
enough to the vertex to make the conical flow assumption valid. 

with the bounderies established and some points in the phys- 
ical plane for which the velocities are known one may vroceed by 
using the Sauer grenhical method for computing strech onic axialiy 
symmetric flow. There are three different geometrical cases for 
wnich slightly different »vrocedures are necesszry. The cases ars: 

(a) from the known solution in tvo points in the free stream 
to calculate the solution in a third point in free stream; 

(bd) from the kmown solution in a point in the free strean to 
calculate the solution in a noint on the body; 

(c) from the lmown golution of a noint on the shock wave end 
a point in the free stream to calculate the solution in a thira 
point on the shock wave. 

“irst, take the case of all points in the free stream. (sse 
Sig. 2). Suppose points P, and - are Inown: to find Ps v4 is 
estimated in the hodogravh plane. Then with meen values, magnitude 
and direction, which are estimated, draw an averaze fas and) se 
line. Their intersection gives F, in the physical nlene from which 
values are obtained for dj ,dh and ~. With these velues one can 
calculate dq and dp from 

se 


dq = Bia on ol f 


dp = an’ & = AN 
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before. The velocities from successive points along the body fall, 
however, in a smooth curve in the hodograph plane making estimation 
easy vith the result that the process generally does not need to 
be repeated. 

Now in obteining a point on the shock wave from ea preceding 
point on the shock wave and a point in the free strean it mist be 
noted that the angle of the shock wave is not known. The shock 
wave is conical only for flow about a cone. For « curved profile 
the shock wave is also a curve which for this method is approxi- 
mated by a series of straight lines, the angles of which must be 
found. To do this the Busemann shock polar is used in the hodo- 
graph plane. The shock poler graphically expresses the relation 
between the velocity direction end megnitude of the downstream 
flow and the angle of the shock wave for a given free stream vel- 
ocity. Referring to Fig. (4) suppose P, and P.. are known; P, is 
to be found, Since P, is on the shock wave it is also on the shock 


polar along with Pi which is known. Zstimete F!' on the shock polar. 


3 
The angle of the ehock wave in the vhysical vlene is the complement 


of the mean of ZOAP, and Z OAP The p line PoP 4s odtained as 


1 3° 
in the free stream case. The intersection of this 0 line with the 
shock wave gives the velues necessary to calculate do. The con- 
struction in the hodograph nlane is made giving the point Ps at the 
intersection of the shock polar with the perpendicular to the } line 


from Ph. Since the shock wave and the ? line fron Po are nearly 


parallel it is often necessary to reveat the procedure. 


HO. 

Then stepping off dq in the 0 a3 direction and dp in the ( ag 
direction perpendiculars are erected as shown in Fig. (2). These 
intersect in a point which will be very close or right on the 
estimated Poe If this is not within the desired accuracy the 
process is repeated until the iteretion shows no change. Table 1 
shows several examples of calculations made in obtaining the values 
necessary to construct Fig. (2). | 

For the case of obtaining a point on the surface of the body 
from a point in the free stream the procedure differs slightly. 
Referring to Fig. (3) suvnose P, is known; Pp is to be found. A 
f line, whose direction is estinated, is drawn from P, until it 
intersects the dDoundary et B. Since the flow is tengent to the 
boundary at all points on the boundary the direction of the flow 
velocity at point B is known, and an estimate of the magnitude of 
the flow velocity is made. From the mean between this estimated 
velocity and the velocity at Pr a Mach angle is calculated. Sime 
{larly a mean value of the angle of flow direction (between Py 
and point B on the surface) is caleuleted. The difference between 
these two angles gives a new ( line intersecting the boundary at 
1P2° The points 7) and ae give the values necessary to calculate 
dq as in the free stream case. As before the construction is made 
in the hodograph plane and the nerpendicular at a distance dq from 
ry along the /) line intersects wo at 1P4- From the mean values of 
Fj and P) draw the corresponding { line fron F, in the physical 


plane. This should intersect the boundary at 1a If not, one 


takes the tangent to the boundary at the new point and proceeds as 








eueT YS TeopsAys 


eS. 


(3) “ANE 


xv ~< -_ RT tS GRE a TEST CR CRC IP TT Pa TORT TINE NIG I TOT BN YES ap ER a I a! ¥? << Ae A 





£2) 





oA’ 7204S 49 aybuv 


4 


aur{.; [LwoysAug eUuBT.; YdBsdoper 


SO9 8 


15. 


Then, Cp = 2 fro foe 
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Z is a function of the strength of the shock wave and thus of the 
shock wave angle, which is known from the Taylor-Maccoll solution 
(the angle of the shock wave is teken at the nose), and of the 
free stream Mach number. 


Given M, and x (see figure below) Z is computed as follows": 
oh v 
“Ma 
Mi, ; 


a a 


ve E | 
———— a / 7 , ey ee 
/ wi J 


/ 


| tea 
aie Ms op (a) 1] 


cotta» Le tee ny Or) 
‘cs 
1 | Ce). 








ten Cre) 4 (re (4) 
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"Durand, ¥. F., Aerodynamic Theory, Volume III, Durand Reprint, 
1943, pp 72388240, 
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CALCULATION OF C,, AND C,, 


After the calculation of the flow velocities 13 completed, 
all of the points of the Mach net bdDetween the shock wave and the 
surface of the body are tabulated, i.e., the geometrical nosition 
of the point and the strerm velocity at the point. (see Teble 11). 
In the calculation of the drag of the nose section, however, only 
the pnoints on the surface of the missile are needed. Given a num 
ber of noints on the surface and the velocities at these points 
the pressure coefficient, *,: is computed for each of these points 
and a plot is made of 0, vs. X. (seo Fig. (8)). From this curve 
velues of ¢, are obtained for regular increments of X and then by 
carrying out the avvrovriate integration (in this case a grephical 
integration was used; see Fig, (10)) the drag coefficient, ©), cen 
be obtained for the nose of the missile. 

The equations used in the actual computations ——r 


below: 


In general Zp Ps. 


end a al. 
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where My is the free stream Nach Number. 
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Across the shock -. however, there is a drop in stagnation 
pressure. This drop is defined by 2 = pz/p, where pz corresponds 
to stagnation conditions behind the shock wave and p, corresponds 


to free stream stagnation conditions. 


16. 
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Ny, the free stream Mach Number is known from the Taylor-Maeccoll 
solution of flow around 2 cone, 4K —___1ocal velocity © . W 


." critical speed of sound ae 


is known at various points on the surface of the body and 


( (4 ‘ a yet) 
5 ~ Sr oe 
Te 
4 
or 
.* 
2 ay 





i = 


< 
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where Mg is the Mach number at any point on the surface of the 
missile. Subdstituting for M, in the equation for p,/p3 

(Y+))- Ve op - ) 
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Substituting these momen iey * ‘pes the expression for c! 
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*Liepmann, H. W., 2nd Pickett, A.E., "Introduction to Aerodynamics 
of a Comressible Fluid", John Wiley & Sons, 1947, p. 26. 
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where the only variable for any one free stream M) is MM", 
The coefficient of drag, Cp, is defined for a three dimensional 


body: } D 


e = oe LL , 2 
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where d is the maximum diameter of the body. 
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= ey? poe 
D 2, f, Cp rn dy 
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This integral is numerically evaluated with ., being known from 


the above eoauation. 


18, 
DISCUSSION OF RESULTS 

The Sauer Method of Characteristics was applied to the nose 
section of the Corporal GZ for free stream Mach numbers of 2 and 4, 
The shape of the missile was given by the parabola r = Is (1 i x ' 
L = 96", The result for a free stream Mach number of 3 was obtained 
from a solution for the same nose section by H. EK. Forster. For the 
case of M = 4 two entirely independent calculations were made to 
check the consistency of results. In one of these solutions the 
nose was assumed to be conical for a distance of six inches from 
the vertex and the constant velocities obtained from the Teylor- 
‘“aceoll solution for e cone were assumed to occur over this vortion 
of the nose, In the second of the commutations the nose was assumed 
to be conicel for e distance of nine inches from the vertex. It 
can be seen from the nlot of pressure coofficients obtained in Fig. 
(€) that the results were consistent. 

The nressure coefficients obtained for M = 2 and M = 4 are 
tabulated in Tables IV, V and VI and ere plotted in Figs. (8) and 
(°), For comparison, values of Coe obtained by the eal Moore 
Method ag amnlied by W. Z, Chien* to Orive Nose Sections, a 
plotted. Figs. (12) and (13). This method makes use of e lineer- 
izing process and a close check of results vould not be expected. 

It is noted however that e closer check is obtained at the higher 
deen numbers. This is consistent with theory. 

Values of the drag coefficients obtained ere: At M = 2, O= 
0.09168, at H= 3, C., = .066, and at = 4, C, = 0.061. These are 


compared to the values obtained from the Karman Moore Method in 


*s3 


te. (24). 


“wave Drag of a Projectile Nose at a Supersonic Velocity by the 
Kerman-Moore Method, ORDCI? Report No. 4~24, Avril 19, 1946. 





2.. 
All of the points computed in the Mach nets are indicated in 
Biers. (5), (6), end (7), and the velocities and shock wave angles 
are tebulated in Tables II and III, (In these tebles the velocity 
ig expressed in terns of the critical velocity, ©", M° = xebocdty) 
The caleul: tion of points end of velocities at these voints 
theoretically could be carried out to any desired degree of accuracy 
denending on the scale used in the graphicel solution. Actually 
the scale used is limited in the hodograph olane by the size of 
drafting eouipment available. In this instance a twenty-four inch 
rule calibrated in 1/50 inch increments (5 units equel 1 inch) was 
used. In the 4 = 4 solution 50 units corresnond to Mf" = 1, and in 
the “ = 2 solution 65 units correspond to M* = 1. ‘These scales 
were chosen so that the full length of the rule would be utilized 
(and thus maxirmm eccuracy) in the hodoegraph plane. In the nhysical 
plane 5/3 units corresnond to 1 inch (f.e. a 1:3 seale). 
Following is an outline of the steps in the solution of 4 = 4 
for the Corporal E nose section: 
l. Taylor-Maccoll solution" 
A, <A value for U/e corresnonding to M = 4 was chosen from the 
preaon, 
B. A tabuler form was set unm and indicated calculations were 
made using increments of g equal to 0.5 cegrees. 
C. The angle of the shock wave was solved for in the Teylor- 
Maccoll region as indicated. 
D0, The next step was to solve for the exact free streen iach 
number. 


"Taylor, G. I. end Maceoll, J. W., “Air Pressure on a Cone Moving at 
high Speeds", Proc. of the Roy. Soc., A, Vol. 159, Fages 276-311. 1933. 
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4. 


8. 


10. 


i. 


nlane in such a position as not to interfere with the Macn net. 


woe 


E. Operations were carried out as indicated in Table VII. 


‘ ty 
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fe 
M™ was converted to the scale in use. 


The nose section wes drawn to the scale chosen. 





A shock polar was constructed for viach number 3.9712 eas indicated 
in Fig. (5). 


Coordinates of the hodogranh plane were drawn above tne physical 


i table of Mach angles and sine squared of the Mach angle was 
made for anticipated valocities. (See Table VIil). 
Y¥-/ he 
—— ae | ae ‘] 
are: - KInAs. A = a tem 
Mm tye? 


r+) 
From Taylor-iaecoll's solution velocities were drafted in the 


hodorfraph rlane and in the vhysical plsne the p lines were 
drawn using average velues of Mach ongle between the. voints. 
(Point 1 was taken ex distance from the nose corresponding to 
6 inches in one case and 9 inches in the other). 

Calculations of the Hach net by the Sauer Ilethod was started. 
it was advisable to caleulete the slone of the nose at the inter- 
section by a f line from the geometry of the body rather than 
to annroximate it graphicelly. 

«nen a point on the shock wave hed Pa be found the shock solar 
wag placed over the hodogranh plane end was removed after the 
noint was located, 

Z wes calculated. (see Table 1%). 

Cy was calculated. (see Table IV), and a plot of C, Vs. x was mad 
“ven increments of X were teken and the corresnonding values of 


r and Gy were tabulated. 


aie, 


12, A plot rO, vs. r, see Fig. (10), was made and the "positive" 
prea ieee this curve was determined. This area was connected 
to the proper scale. (In this case the conversion factor 
was (2) X 0.4 ). 

13 The value for Cp = fe x Aree was calculated. 

It was found thet the most efficient method of conatructing 
the kach net was to work in pairs. One verson drafted while the 
other verforned the calculations. Fron three to ten points per 
hour were cormleted after nroficiency in the mechanical processes 


was obtained, The average working sneed was about five points ver 


hour, 
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In Column 5, 
50 units corres»ond TABLE I 
to M® =] 


SAMPLE TABSLI SHOVING DETAILAD 
CALCULATIONS BOR MACH NRF 

















a. a ve Jf ix —an< 
S&S « a = 

189 102.05 9 15 16.35 

169 = 68 101.87 9 22 17 9 0870 16.58 

169 - 61 102.05 9 39 17 4 0862 17.22 

5S ta ) . 7 5 | 
170 102.46 8 60 15.87 

170 = 69 102.25 9 2 16 57 0849 16.11 

179 = 62 102.26 9 20 16 56 0849 16.79 

70 102.46 8 5,87 
171 103.12 8 15 14,84 

171 = 70 102.79 8 32 16 38 0819 15.35 

171 = 63 103,22 8 60 16 24 0797 15.96 

7} 103.16 8 25 15.10 
172 103.88 7 §8 14.46 

ee = 71 103.57 s a 16 04 0766 14.78 

; 1 7 14,44 

173 101.54 8 44 15.48 

17> = 65 101.37 9 14 17 37 0899 16.10 

73 101.54 8 44 15.48 
74 65 = 73 101.37 9 =? 16.70 
175 101.48 8 55 15.82 

175 = 74 101,42 -9 #2 7? 27 0899 15.96 

175 = 66 101.40 9 27 1? 37 0899 16.32 

5 101.50 9 6 16.00 
176 101.72 8 52 15.80 

176 = 75 101.61 8 58 > 20 0887 15.90 

176 = 67 101.69 9 8 17 16 0881 16.23 

76 7 101.72 8 52 15.80 
177 101.42 9 0 16.04 

177 - 76 101.57 8 56 17 2l 0889 15.92 

177 = 68 101.56 9 15 17 2 0889 16.43 

77 101.73 8 50 15.90 
178 102.00 8 48 15.61 

176 = 77 101.87 8 49 17 10 0872 15.76 

178 = 69 102.03 9 O02 17 O05 0864 15.98 

76 102.00 8 48 15.61 





Point 73 is a point on the shock wave. Point 72 is a point on the missile. 
All other points are free stream points. P denotes the column for estimated 
points. The individual estimated points are designated by ,( ), o( ), etc. 
i$ ) - ( ) denotes average values between the indicated points. 


1. a oe a 


TABLE I 


—_ — + 


(Oont 'a,. ) 


53. 


eo 


air 











22440 

22.265 Sewe ° c6 «61 -7 47 

20.87 485 26 «43 -? 25 
69 22440 

1.855 

e2eLl 3.96 244 29 59 -7 35 

20.31 6.90 4835 26 «16 -7 36 
70 21.83 

20.90 

21.36 6.62 427 ap 10 <8 6 

19.45 6.82 445 25 14 -7 $4 

20,90 

19.96 

20.43 6,90 381 24 13 -7 53 
fe 19.96 

50.354 

27.18 14.0 » 745 26 41 Se Ts 
73 30.34 
74 27.18 

26.82 

26.95 2.38 127 26 29 -8 25 

ar 6,92 40 £6 54 8 O 
79 26.82 

26.40 

26.61 2.83 015 26 -8 22 

24.89 6.98 » 404 26 6 -8 8 
76 26.40 

26.00 

20420 2.65 143 26 17 -8 25 

24,45 6.90 e412 ce 86 -& 6 
7 26,00 

25.50 

25.75 3.75 . 200 ce oo -8 al 

23.95 7.00 405 a5 «2@ -8 3 


78 29290 


Paaly bad (2) 
SLOT ¥RLOCITIES FOR 


In Colum 2 % = 5.9712 
566 units corresvonds San i 
to "7 = 1 
.oint - a Foint Bm 
ge 2 
1 98, } 
; =| a 46 nh 
3 98.81 14 28 as ae 
4 93.05 1S A i2Z = an 
5 99,40 1 [Gy = ee 
6 99.62 i a 0 ae or 
? 99.92 15 34 s ie 
8 99.92 13) 83 - pers 
9 99.92 15 0? os ee 
10 99.95 1? ~—s 88 a ates 
11 99.98 18 387 = Sore 
12 99.98 14 37 Ze ee 
13 100.17 ll 29 a a 
14 100,02 2 6 6 Ee ead 
15 99,96 15 ae a EeaL 
16 100.08 14 02 Be 10) ee 
17 «100.38 10 41 ne oma 
12 100.18 ll 25 Ee oe 
19 100.10 © 12 £490S o aes 
20 100,22 13 00 of ete 
oT 100.53 13 25 Ne ae 
22 100.00 1l 805 = eng 
a5 100,48 10 15 2? oe 
24 100.24 10 3640 ae ae 
25 100.60 10 «235 ha ap 
26 100.45 10 53 gy ee 
27 100.36 li 2 SS ae 
28 100.45 12 sh Oe 
2g 100,74 i & iy ee 
30 101.17 iz (« he ae 
31 100.82 9 48 a ea 
= 100.61 10 2 a ae 
oe 100.56 11 O - a 
34 100.55 mM | 640 i oe 
35 100.79 i 60S ~ 1 Oats a 
36 101.11 i «i Y peta 
SF 101.25 ~~ YF = ooo ae 
33 19] 39 9 10 We ioe 
39 101.06 9 29 = Se 
aD 100.79 10 0 oe en 
4} 100.68 10 34 we wore 
4° 100.69 ie iE ee: 
13 100,88 11 «413 - Sore 
Ad. 1Oieat 7 ll @ se ttc-ae 
45 101.37 1l 365 - ee 


90 104.61 





“i” 


105.24 
101.96 
101,83 
102.00 
102.15 
102.31 
102.70 
103.35 
103.85 
104.65 
105,30 
106.44 
102.13 
102.15 
102.43 
102.78 
103.41 
103.98 
104.76 
105.28 
106.33 
107.49 
102.45 
102.29 
102.44 
102.56 
102.93 
103.43 
105.98 
104,64 
105.°5 
106,28 
107.25 
107 as 
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In Column °, 
59 units corresnonds 356. 


to #* = 1 
TABLE Ii (h) 
FLOW VELOCITIES FOR M # 3.9712 
Sor II 
( ovat ae 7 Point M* 
oy = 
98.60 17 +08 46 101.10 
2 98.75 15 18 4? 101.45 
3 98,97 13 ««42 48 101.60 
4 99.49 2 Sey 49 102,00 
5 99,50 16 05 50 101.00 
6 99.55 3 85 51 101.00 
7 99,70 “6 52 101.10 
8 99.90 le 10 53 101,40 
9 99,40 14 00 54 ) 101.55 
10 99,80 15 2 55 101.85 
11 99,62 1? 39 56 102.00 
12 100,10 ll 35 57 101.25 
13 99.55 12 60 58 101.12 
14 99,85 13. «50 59 101.10 
15 109,00 14 49 60 101.25 
16 99,86 ll 6-28 61 101.50 
1? 100,15 10 30 62 101.60 
18 99.85 li 45 as 101.90 
19 100,00 12 «(15 64 102.00 
20 100,00 137 65 102.60 
om 100.15 iy Te 66 101.30 
90 100,05 10 55 67 101.40 
- J 03 100,35 10 15 68 101.65 
24 100.10 ll 9625 69 101.70 
23 106.20 1l 8645 70 102.10 
26 100.25 12 «(3s 71 102.20 
2? 100.35 13 28 72 102.60 
28 109.90 isle 73 102,00 
29 109.19 1. «Ge 74, 101.60 
zO) 109,35 190 15 75 101.45 
3] 100,35 12 8650 76 101.60 
eo 100,40 11 & 27 101.70 
23 100.50 12 #15 78 101.90 
34 100.90 12 79 102.20 
35 101220 13 4&2 80 102.350 
36 100,90 30 82 102.65 
37 100,53 28 82 103,00 
38 100.60 45 83 103.60 
39 100.60 35 B4 101.80 
40 101,00 45 85 101.95 
4] 101.40 0 86 102.00 
42 101.59 25 87 102.25 
AB 101.72 56 88 102.40 
a4 100,30 ie 39 162 70 
45 100,90 55 90 103,05 





Point 


1i4 
de 
26 
LL? 
Tee 
stats 
1°06 
17) 
je 
eS 
1?4 
ves 


139 


Mm" 


103.60 
104.15 
102.10 
101.95 
Mee 360 
102,05 
102.30 
102.95 
103.95 
104, 00 
104,80 
L0?.10 


102,20. 


102.00 
103,00 
103.55 
104,05 
104,85 
105.55 
102.50 
LOZ a0 
102. 40 
102.60 
103.10 
103.65 
104,10 
104, 80 
105.50 
106.30 
107.35 
L0?.65 
103,20 
103,70 
104.16 
104,95 
105.45 
LOG. 30 
107,30 
102,60 
102.48 
102,85 
103.30 
103,80 
104.35 
105,00 
195.65 
106.40 
107.40 
tO. 00 
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TABLE II (bd) Cont'd. 


J 


In Colunn a 38, 
65 units corresxonds 


to “M*= 1 
TABLE Ili 
FLOY VSLOCIPIsS FOR H = 2,0039 
Point hm J Point i{* 

2 - Q 

1 94.90 13 46 46 99.38 6 
S 95.65 ll 19 47 99.60 6 
3 96.38 9 7 48 99.05 t 
mae 97.28 8 a 49 oF ..70 a 
5 98,01 7 28 50 98.70 8 
6 99,44 6 &5 ol 98,65 g 
7 98.94 6 28 92 98.50 16 
& 99.19 6 Os 53 99.00 ll 
9 96.05 16 0% 54 99.69 Le 
10 26.200 i2 40 55 99.58 c 
11 95.20 15 4s 56 99.95 D 
2 fO.05 LO 35 oT 99.40 6 
uS 96.20 12 35 58 99.10 ? 
re 95,35 15 30 59 96, 20 7 
15 97.40 8 55 60 oe 8 
16 97,00 10 ic 61 99.10 9 
1? 96.30 11 58 62 98 WLS 9 
18 96.50 1S OF 63 99.70 Le 
19 97,30 8 OC 64 1GOe Zo 11 
20 97.05 8 45 e 100.25 5 
ol 97600 10 05 66 99.80 6 
ye 96.385 le Q1 67 99.50 G 
25 97.30 14 28 68 99.55 ? 
24 sie yore 7 25 69 99.45 ? 
25 98.10 8 00 70 99.45 8 
26 97.50 9 05 71 9950 8 
2? 97.20 10 35 Ve 99.90 9 
28 97,40 12 20 73 100.20 10 
29 97,55 14 05 74, 100.8 . oe 
30 oa. 7 6 50 75 100,80 “ 
31 98.65 ? le 76 100, 43 5 
a 98.15 8 00 77 100.10 5 
30 97.80 9 15 78 99.80 6 
34 97.95 10 20 79 99.50 0 
35 98,00 li 30 80 99,40 7 
36 98.50 13 28 81 99.85 8 
37 99.10 6 30 8> 100.138 8 
38 99,00 6 50 83 101.00 9 
ag 98.50 7 40 84 101.00 10 
AO 97.95 8 40 85 100.62 4 
41 98.15 9 40 86 100.30 5 
42 98.20 ZO 25 87 100.10 5 
43 98,50 le OO 88 99.80 6 
44 98,60 13 10 89 99.75 ? 
45 99.78 5 40 90 100,05 7 





114 
115 


117 
118 


120 
lal 
Lae 
123 
104 
1Zo 
126 
127 
128 
128 
130 
131 
132 
133 
134 
135 
136 
137 
138 
139 


M* 


100,30 
100.95 
100.30 
101.80 
100.50 
100.30 
100,20 

99,85 
100.10 
100,30 
101.20 
100.95 
101.55 
102.20 
100.75 
100.62 
100.50 
100.35 
100.30 
100,40 
100.55 
101.35 
101.20 
101.65 
102.20 
103.05 
100,62 
100,45 
100.45 
100,50 
100.60 
101.40 
101.35 
101.65 
102.00 
102,80 
103.75 
101.25 
100.83 
100.75 
100.80 
101.40 
102.05 
103.60 
104, 80 
100,04 
100.70 
100.80 
101.35 
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TABLE III (Cont'd.) 


J 


Point 


140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
151 
152 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
168 
169 
170 
171 
172 
173 


102.00 
103.45 
104,45 
105.65 
101.20 
101.35 
101.80 
102.35 
163.75 
104, 80 
105.90 
107. 80 
101.24 
101.65 
101.44 
101.35 
101.40 
101.80 
102.35 
103.70 
104.65 
105.70 
107.45 
108.20 
101.65 
101.75 
102.10 
102.60 
103, 80 
104.70 
105.70 
107,40 
108.25 
109.00 


to 
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roint 


bes 
16 
ol 


a” 
46 
54: 
633 
72 
82 


102 
Deo 


L2¢ 


FOR CALCULATION 
Cc FOR HM = 3.9712 
(élite TO FIG. (6). 


Coordinates 

ie ~ 

lo 4.56 

18,83 5.48 

COede 6.45 

26.285 7.49 

52.62 9s 

37.96 10,45 

43.19 11.67 

48.74 12 seo 

55.48 14,34 

63.05 LO. fe 

75.20 18.00 

&8.17 19.96 

101.97 21.75 
116 2¢3 2312 
136.40 24,46 
LT. oe 24.98 
182.10 25.00 


PALLY IV 
DAPA SHEET 


50 * 


93,60 
99.93 
99.81 

99.98 
100.08 
100.53 
101.17 
Lola 
LOT @ 72 
102,50 
103,35 
103.75 
104.60 
105.24 
106.44 
107,49 
107.18 


if 
if® 


1.9720 
1.9086 
1, 982 
1.9996 
229016 
2.0106 
2.02354 
2.0250 
2.0358 
2.0500 
9, 0666 
2.0750 
2.0280 
221048 

2.1288 


> ee 


221436 


Cy 
2209 
~1621 
01634 
» 1604 
01596 
1419, 
lela 


Lee. 


1032 
0833 
0619 
0519 
0330 
0201 
-.0016 
=-.0178 
-.0133 


Fat 

he 
~0.66500 
@ O.9, ef 


40, 
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TABLE V 
DATA SEETT FOR CALCULATION OF oe FOR if = 3.9712 


CORRESPONDING TO FIG. (5) 


x R 50 M® M® G ro 
P P 
10 3.03 98.6 1.972 »28O9 «6550 
12760 3,80 g9 6 1.980 ©1957 « 7436 
16.20 4.75 99.70 1.994 1701 8079 
e0.6p ‘6,00 99.80 1.996 - 1666 9996 
24.85 7,10 100.00 2.000 ~1597 1.1333 
30.20 68.45 100.15 22003 1546 1.3063 
bD.-60 8.75 100.90 2.018 ©1301 1.2684 
41.15 11.10 101.20 7 «020 1207 1.33597 
47.20 12.50 101.50 2-050 1116 1.5950 
52.80 13.80 102.00 2 040 0971 1.5400 
58.30 14.95 102.00 2 eQ40 0971 1.4516 
65.90 16.40 102.60 22052 - 0806 1.d218 
75.70 17.70 103.00 2e060 «0702 1.2425 
83.30 19.30 103.60 eum 0554 1.0682 
93.00 20.65 104.15 2.083 . 0428 . 8338 
107.80 22.30 104.80 2.096 0287 - 6400 
122.80 23.70 105.55 Ge Lae ~Q132 03128 
139.10 24.60 106.30 22126 0007 0172 
157.65 24.85 107,30 2.146 -.0151 -0.5752 
180.80 25.00 107.50 22150 -.0179 -0.4475 


TABLE VI 
DATA SHE? FOR CALCULATION OF C,, FOR MN = 2.00389 


CORRESPONDING TO FIG. (9) 


Point x r 65 M* M® C., 

0 10.00 3.08 94,57 1.454922 eo8239l 

9 12.40 3-80 96.05 1.477692 24400 
11 15.00 4.50 95.90 1.475384 004292 
14 18.00 5280 95.85 1.474615 ~20051 
18 22,80 6.65 96.50 1.484615 o2d419 
23 27.90 7.95 97.30 1.496923 021442 
29 31.50 8.85 97.55 1.500769 20831 
36 37.60 10.30 98.50 1.515538 218513 
44 40.50 11.00 98.60 1.616923 ~18299 
54 47.30 12.60 99.65 1.533076 15822 
64 53.80 14.00 100,35 1.543846 14205 
74 60.50 15.40 100,80 1.550769 ~13178 
84. 67.00 16.56 101.00 1.553846 12724 
94 74,00 17.80 101.80 1.566153 ~ 10934 
104 80.80 18.90 102.20 1.572307 o20053 
116 88.70 20.10 103.05 1.685384 08205 
127 96.40 21.10 103.75 1.596153 06732 
1356 110,00 22.60 104, 80 1.612307 .04528 
143 123.60 23.90 105.65 1.625384 02805 
151 139.10 24.70 107.80 1.658461 -+,.01375 
163 153.20 25.00 108.20 1.664615 «,02122 
173 168.80 25.00 109,00 1.676922 -,03592 





TABLE VII 
TAYLOR-MACCOLL: KF = 2.0039 


VELOCITIES AT A MISSILZ TIY FROM PAYLOR-MACCCLL SOLUTION 
D =v u M* XO tlt -@ 
; ae 6 lO 

17,08 0 ~89700 1.4549 3712 0 0 
1, 00 03852 .59635 1.4564 37 8 .06459 3 42 
21.00 ~07538 .59435 1.4600 $6 56 .1°7683 714 
33,00 10992 .59111 1.4651 36 44 .18595 10 32 
20.0 ~14283 .58670 1.4715 36 26 .24344 13 41 
76.5 ~16674 .58265 1.4768 3615 .28618 15 58 
28.0 ~L9017 .57798 1.4829 35 56 .32902 18 13 
29.9 2lde8 .57210 1.4892 35 39 .37241 20 26 
51.0 ~23626 6.56682 1.4966 35 20 .41682 22 38 
32.0 ~-26162 .56257 1.5019 385 8 ,.44727 24 6 
33.0 226713 .55804 1.5078 34 51 .47869 25 35 
34.0 28795 .55325 1.5144 34 $4 .51145 27 5 
35.0 29938 .54818 1.5222 6414 .54613 28 38 
36.0 -S1711 .54280 1.55319 33 50 .58421 30 18 
35.333, -30596 .54639 1.5260 3405 .55997 29 15 


p = geonetrical angle selected for Taylor-"accoll solution 


p = arctan 


anrle of shock wave 


Mh ee elo 
oo A 


pas 
- % 


J 
Ooaonnenwvd ow OO 


5 
18 
46 
<8 
19 
o2 
47 
ae 
20 
55 


“2 


QO5 





GABLE VIII 


VELOCITIES 4f MISSILS TIP FROM TAYLOR-MACCOLL 
SQLUPION M = 3.9712 


: f 
Dp y” ae ue “vw a p J, 


= oe 


Ces oO 4 

18 00066 .65448 1.972953 1904 .031656 1 49 16 
We ,006117 .651540 1.97637 18 67 .096893 5 32 14 

a .01656 .64571 1.98324 1844 .16013 9 OS 12 
63.8 .02784 .65950 1.99081 1627 .s0868 Wm 47 21 
raealo 
19 .00427 .65256 1.97507 1859 .06475 3 42 15 
2) .0L072 .64898 1.97958 1851 .12853 7 19 ig 
mS ,0371 .64175 1998881 i986 .iofe, 10 Ga Ge 


p = geometrical angle selected for Taylor-Maccoll solution 


/ 
p = arctan 
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